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We consider a general class of (intersecting) loop models in d dimensions, including 
^S| ! those related to high-temperature expansions of well-known spin models. We find 

p^ ' that the loop models exhibit some interesting features - often in the "unphysical" 

region of parameter space where all connection with the original spin Hamiltonian is 
C\J , apparently lost. For a particular n = 2, d = 2 model, we establish the existence of a 

^^ I phase transition, possibly associated with divergent loops. However, for n ^ 1 and 

Qs ' arbitrary d there is no phase transition marked by the appearance of large loops. 

0^ . Furthermore, at least for d = 2 (and n large) we find a phase transition characterised 

-y ' by broken translational symmetry. 
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d ' 1 Introduction 



In the recent years there has been much interest in various loop models. Loop 
models are graphical models defined by drawing closed loops along the bonds 
of the underlying lattice. The loops may come in n different flavours (colours). 
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No two loops can share a bond, while sharing a vertex is generally allowed. 
Explicitly, the bond configurations are such that each vertex houses an even 
number - possibly zero - of bonds of each colour. Each loop configuration is 
assigned a "weight" that depends on the number of participating vertices of 
each type. In the cases of interest these weights are actually positive hence, 
at least in finite volume, they define a probability measure on the set of all 
loop configurations. Thus, for a finite lattice the loop partition function may 
be written as: 

Z = ^ i^VrVa"^ . . I/y ^ , (1.1) 

g 

with the sum running over all allowed loop configurations Q. Here b is the total 
number of participating bonds, rrii [i = 1, ... ,V) is the number of vertices of 
type i and Ui is the corresponding vertex factor.[^ This definition is slightly 
different from the one typically found in literature (c/. Refs. [1,2]) since it 
also includes the bond fugacity R. Although strictly speaking it is not needed 
(since the bond fugacity can always be incorporated into the vertex factors), 
we find it convenient to keep i? as a separate parameter. We remark that 
by relabeling the empty bonds as an additional colour, these models may be 
formally regarded as "fully packed" . 

The reason loop models have been extensively studied is because they appear 
quite naturally as representations (often approximate) of various statistical - 
mechanical models. These include, among others, the Ising model (this ap- 
proach dates back to Kramers and Wannier [3] and was later used to solve the 
model exactly [4,5]), the Potts model (polygon expansion [6]), 0(n) spin mod- 
els [7-11], 1-D quantum spin models [12], a supersymmetric spin chain [13], 
the g-co louring problem [14,15] and polymer models [16,17]. 

Here we consider the loop models explicitly related to the high-temperature 
expansions of the standard 0(n), corner-cubic (AKA diagonal-cubic) and face- 
cubic spin models. This is, in fact, the same set of models that was treated in 
Ref . [7] . However, in this paper, we provide a careful treatment of the large 
n cases - and we treat the standard (i-dimensional lattices. As a result, we 
arrive at quite unexpected results concerning the behaviour of these models 
in the high fugacity region. 

In particular, despite the considerable attention the subject has received, most 
authors (with certain exceptions, e.g. [3-5,9,13]) chose to consider models 



^ Many authors consider an additional factor of the form F^F2 ■ . . F^ where Fi 
is a "loop fugacity" and /j is the number of loops of the i-th colour. Although the 
objects li are unambiguous when self-intersections are forbidden, in the general case 
they are not easily defined. Nevertheless, the essence of such a term - at least in 
the case of integer Fi - is captured by the introduction of additional colours. 



where only loops of different colours are allowed to cross each other (if at 
all). On the other hand spin systems (in the high-temperature approxima- 
tion) naturally generate self-intersecting loops. In order to avoid this issue, 
an exorbitant amount of work has been done on lattices with coordination 
number z = 3 (e.g. the honeycomb lattice), where loop intersections simply 
cannot occur. Overall this approach appears to be justified since one is usu- 
ally interested in the critical properties of the underlying spin systems. Indeed, 
consider the archetypal n-component spin system with |Sj|= 1 and let us write 
exp(AX](jj) Sj- Sjj ~ Y\{i,j) (1 + -^Sj- Sj). Although as a spin system the right 
hand side makes strict sense only if |A|< 1 (the "physical regime"), the as- 
sociated loop model turns out to be well defined for all A. Since the systems 
can be identified for |A|^ 1 it can be argued that the critical properties of the 
spin system and those of the loop model are the same and are independent of 
the underlying lattice. 

Notwithstanding, for n ^ 1 any phase transition in the actual spin system is 
not anticipated until temperatures of order 1/n (i.e. A ~ n), which we note is 
well outside the physical regime of the loop model. At first glance this appears 
to be borne out: the natural parameter in the loop model (as well as in the spin 
system) seems to be \/n. Thus, the loop model could, in principle, capture 
the essential features of the spin system up to - and including - the critical 
point. 

We have found such a picture to be overly optimistic. Indeed, depending on 
the specific details, e.g. the lattice structure, there may be a phase transition 
in the region 1 ^ |^|^ ^ (specifically, A ~ n^/^), well outside the physical 
regime but well before the validity of the approximation was supposed to 
break down. Furthermore, it would seem that both the temperature scale and 
the nature of the transition (not to mention the existence of the transition) 
depend on such details. Finally, we shall demonstrate that in contrast to their 
spin system counterparts, the large-n models have no phase transition - for 
any value of bond fugacity - associated with the formation of large loops (i.e. 
divergent loop correlations). 

The structure of this paper is as follows. Section 2 is dedicated to the de- 
scription of the spin models and their connection to the loop models. Specific 
results for those models with the two-dimensional spin variable (n = 2) are 
presented in Section 3. Finally, Section 4 contains the discussion of reflection 
positivity as well as some results concerning phase transitions in the large n 
case. 



2 n-component models 



2.1 Oin) model 



Let us start by considering the 0(n) model on some finite lattice K d 7/ 
defined by the following partition function: 

Z = Trn(l + AS, -S,) (2.1) 

with Sj G M"^, |Sj|= 1 and Tr denoting normalised summation (integration) 
over all possible spin configurations. The corresponding loop model is readily 
obtained along the lines of a typical "high-temperature" expansion. We write 
Sj -Sj = Si Sj + . . . + Si S^ and define n different colours (each associated 
with a coordinate direction of the O(n)-spins). Expanding the product, we 
have n choices for each bond plus a possibility of a vacant bond. Thus, various 
terms are represented by n-coloured bond configurations: Q = {Qi, . . . ,Qn) 
with Oi denoting those bonds where the term S^ Sj has been selected. 
Clearly, the various QiS are pairwise (bond) disjoint. Thus, for each Q we 
obtain the weight 

Wg = TT n ASf^S-f^.. n ASf^S-f^ (2.2) 

{i,j)&Gi {i,j)&Gn 

On the basis of elementary symmetry considerations it is clear that Wg 7^ if 
and only if each vertex houses an even number (which could be zero) of bonds 
of each colour. Once this constraint is satisfied, we get an overall factor of A^*^^) 
- with b{Q) being the total number of participating bonds - times the product 
of the vertex factors obtained by performing the appropriate 0(n) integrals. 
The details and results of these calculations are presented in Appendix A. 
In general, it is seen that the vertex factors depend only on the number of 
participating colours and the number of bonds of each colour emanating from 
a given vertex, i.e. not on the particular colours that were involved nor on the 
directions of these bonds. 

In the case of a square lattice {d = 2) we have only three main types of (non- 
empty) vertices: those where two bonds of the same colour join together, those 
with two pairs of bonds of two different colours and those with four bonds 
of the same colour. These have weights of 1/n, l/n{n + 2) and 3/n(n + 2) 
correspondingly. Rescaling the bond fugacity from R = X to R = X/n we 
arrive at the vertex weights ui = 1, z/2 = n/{n + 2) and z/3 = 3n/{n + 2). 

The factor of 3 relating 1/2 to u^ has an interesting interpretation (which, as 
shown in Appendix A, turns out to be quite general). Indeed, each vertex of 
the third type may be decomposed into three different vertices as shown in 



Fig. 1. Each of the new vertices is now assigned equal weight, which is also that 
of V2- We thus split each Q into 3™"^ different graphs - each of equal weight - in 
which every vertex with four bonds now provides explicit instructions relating 
outgoing and incoming directions of an individual walk. 



+ 



+ 



Fig. 1. Decomposition of a type 3 vertex into three new vertices in the 
two-dimensional 0(n) loop model. 



Hence in every such graph (now defined with the walking instructions encoded 
at every vertex) the individual loops are well defined. Furthermore, changing 
the colour of any loop does not change the weight of the graph. Thus we may 
write 






(2.3) 



where the summation now takes place over all configurations /C of colourless 
loop graphs in which every vertex housing four bonds is resolved by "walking 
instructions" , and i is the number of such loops (being now defined completely 
unambiguously). In addition to the advantages of a manifestly colourless ex- 
pression, the above permits continuation to non-integer n. 

We conclude this subsection with the following series of remarks and observa- 
tions. 



As shown in Appendix A, such vertex decomposition works, in fact, for 
an arbitrary lattice in an arbitrary number of spatial dimensions (with the 
proper weights for vertices housing 6, 8, etc. bonds). 

Notice that only |A|< 1 region of the parameter space is "physical" (in a 
sense of the underlying Hamiltonian: —PH = Y.{i,j) ln[l + ASi-Sj]), while for 
|A|> 1 one presumes, no spin Hamiltonian can be written at all. The corre- 
sponding loop model, however, makes perfect sense in the entire parameter 
space. 

If we consider a 2D XY model (n = 2), we notice that the factor 2^^^ ^ in 
Eq. (2.3) can be obtained by assigning directions to the colourless loops. The 
above decomposition of type 3 vertices makes this procedure unambiguous. 
Having done that, we can turn this model into a random surface model by 
assigning heights to the plaquettes in such a way that a plaquette to the 
right of a directed bond is always one step higher than the plaquette to the 
left. Not surprisingly, this random surface model turns out to be identical to 



the one obtained by the standard means of Fourier-transforming the original 
weights in Eq. (2.1) (c/. Refs. [18,19]). 

Finally, it is worth mentioning that in the fully packed limit R ^ oo, the 
n = 2 loop model on the square lattice (with the colour degrees of freedom 
being replaced by assigning directions to the loops) turns out to be nothing 
but the square ice model (i.e. the six-vertex model with all six weights being 
equal - see Ref. [6] for a definition of this model). The mapping between the 
vertices of these models is shown in Fig. 2. We remark that the perspective 
of the ice model (and, for that matter, other six-vertex models) as a two 
colour loop model provides additional flexibility in the analysis of these 
systems. These issues will be pursued in a future publication. 
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Fig. 2. Mapping of a fully packed 0(2) model onto the square ice modeL 



2.2 Corner-cubic model 



We now consider the following "discretised" modification of the above 0(n) 
model (given by Eq. (2.1)): 



Trn 



{i,j) 



X 



n 



(1) (1) , (2) (2) , (n) (n) 

(J- V] + 0-- 'ay ... + al 'a) 



[2.A) 



with ai = ±1. For small values of A this model may be viewed as a high- 
temperature limit of a corner-cubic model. Indeed, it describes an interaction 
of the type in Eq. (2.1) where spins Sj are allowed to point at the corners of 
an n-dimensional hypercube (with the origin being placed at the centre of the 
cube). 



Mapping it onto an n-colour loop model is almost identical to the 0(n) case. 
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with the only difference being the vertex factor: (fj ) • • • {(^l"' ) = 1- We 
can choose to associate the weight of i? = \/n with each bond, thus making 
all vertex weights z/j to be equal to unity. In other words, the resulting loops 



in this model do not interact with each other via vertices (there is still a hard 
core bond repulsion, however). The partition function is then simply 
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Fig. 3. A fragment of a possible two-dimensional loop configuration for an inter- 
secting loop model. Three different possible colourings (out of the total of 32 in the 
case of n = 2) for a loop model of corner-cubic type are represented here as (a), (b) 
and (c). For the face-cubic type, only monotonous clusters like (a) remain allowed. 



2.3 Face-cubic model 



Finally, let us examine a different model with cubic symmetry given by the 
following partition function: 



("),,(") 



Tr n [l + A {^4'^y^^ + v?^v^ . . . + v^v^ 

{iJ) 



(2.6) 
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Here u\'^' = 0, ±1, and for a given site i exactly one of u\'^' {k = 1, 2, ... n) 
has a non-zero value. In fact, one may think of m's as components of an n- 
dimensional unit vector that is only allowed to point along the coordinate axes 
(or from the centre to the faces of an n-dimensional hypercube - thus the name 
face-cubic). While the corner-cubic model described earlier had 2" degrees of 
freedom per site (the number of corners of a hypercube), the present model 
has only 2n such degrees of freedom (the number of faces). 

Once again, the corresponding loop model is obtained by performing multi- 
plication in Eq. (2.6) and then summing the resulting terms over all possible 
values of -u's. But since for each site i only one of the spin components -Uj 7^ Q 
at a time, we notice that no terms that mix different /c's are allowed. In terms 
of resulting loops this means hard-core repulsion of different colours: only loops 
of the same colour can share a vertex. The vertex factors are now: zero for 



any vertex with multiple colours and 1/n for vertices with two or more bonds 
of the same colour; the bond fugacity is given by -R = A. 



3 Results for the n = 2 case 



3.1 The n = 2 models with cubic symmetry, Ashkin-Teller and random sur- 
face models 



In this section we shall restrict our attention to the models with cubic symme- 
tries. Firstly, let us slightly change the notations for convenience: let al = ai 
and al = Ti for the corner-cubic, while u\ = Ui and u\ = Vi for the 
face-cubic model. The corresponding partition functions are then written as 



^cc = Tr n 

{id) 



1 + - [aiaj + TiTj^ 



(3.1) 



and 



Zfc = Tr n [1 + a (uiUj + v,Vj)] . (3.2) 

While we have seen that the loop models generated by these partition functions 
are very different, the spin models themselves turn out to be identical. Indeed, 
Eq. (3.2) is obtained from Eq. (3.1) by the following transformation: Ui = 
((Tj + Tj)/2, Vi = {ai — Tj)/2. This is equivalent to a 45° rotation in the spin 
space (along with a v2/2 rescaling) and is very specific to the n = 2 case. 
In turn, both models are equivalent to the Ashkin-Teller model [20] with a 
particular choice of parameters that will be detailed in Section 3.3. 

The two-colour loop models generated by Eq. (3.1) and Eq. (3.2) are given by 
the following sets of parameters in Eq. (1.1): R = A/2, with all vertices having 
weight one, and R = X, with all multi-colour vertices given weight zero and 
all other non-empty vertices given weight one-half respectively. 

Turning our attention to the particular case of two spatial dimensions, we 
remark that in the former model one can sum over all possible colourings to 
obtain the following result for the partition function: 



■^^KG') 



2^(^') (3.3) 



with b{Q') being the total number of occupied (colourless) bonds and f{G') 
being the total number of faces in the clusters they form. The number of 



faces is the minimum number of bonds that one must remove in order for the 
remaining clusters to be tree-hke. For example, the cluster in Fig. 3 has five 
faces, while it can at most consist of four loops. Curiously enough, this result 
appears to have no simple generalisation for n > 2. 

It appears that the two-dimensional loop model derived from the corner-cubic 
model can not be mapped directly onto a random surface model. However, 
the other loop representation (the one obtained via expansion of the face- 
cubic model) does correspond to a random surface model. Indeed, consider the 
following "recipe": take a loop configuration generated by a particular term 
in the expansion of Eq. (3.2). Let red be the colour of loops originating from 
m's, while blue corresponds to f 's. Take all plaquettes at the outermost region 
to be at height zero (these plaquettes are said to form a substrate). On this 
substrate we have clusters of loops. The outermost boundaries of these clusters 
are themselves closed loops. The plaquettes immediately adjacent to these 
boundaries are assigned the hight of +1 if the loop forming a boundary is red, 
or —1 if it is blue. These plaquettes, along with any other plaquette accessible 
from them without crossing coloured bonds are said to form a plateau and thus 
have the same height .[^ Inside such a plateau region there may be other loop 
clusters, that may or may not touch the boundary of a plateau (only corners 
are allowed to touch, since no bond sharing between the loops is possible). 
Every such cluster is now treated in the same way: its boundary defines the 
"secondary" plateau with the height being that of the "primary" plateau ±1 
depending on the colour of the boundary. This procedure is repeated until all 
plaquettes are assigned their heights. As an example, consider the cluster in 
Fig. 3(a), which may now only consist of the bonds of a single colour. If this 
were a red cluster, the heights would be +1 for the plaquettes 1-4 and +2 for 
the plaquette 5. 

In fact, this description is essentially identical to that given in Ref. [21] in 
the context of wetting transition with the only difference that we allow for 
two-coloured clusters instead of single-coloured, and therefore the heights in 
our case may be both positive and negative.^ 



^ The strict definition is as follows: plaquettes A and B are said to belong to the 
same plateau if and only if there exits an unbroken path along the bonds of a dual 
lattice that connects the centre of plaquette A to the centre of plaquette B without 
crossing a single coloured bond of the direct lattice. 

^ This random surface model, however, has a few important differences with those 
of a more conventional type (like the one obtained for the 0(2) case). Firstly, due to 
an irreducible four-leg vertex factor, it cannot be described by a nearest-neighbour 
Hamiltonian (i.e. a Hamiltonian that depends only on the height difference of neigh- 
bouring plaquettes). Secondly, since no directions are assigned to the loops sepa- 
rating the plateaux, there is no way of deciding on the sign of their relative height 
difference without going through the necessary construction steps starting from the 
outside. By contrast, the 0(2)-related random surface model can be constructed 



The important feature of this random surface model is that it must have a 
phase transition whenever the underlying Ashkin-Teller model undergoes a 
transition. 



3.2 Random cluster representation 



Let us derive yet another graphical representation for the n = 2 (Ashkin- 
Teller) model considered in the previous section, this time it will be a random 
cluster representation closely resembling the FK representation for the Potts 
model. We start from Eq. (2.4) (with n = 2), and with the help of the identity 
(TjCTj = 2(5o-jo- — 1 rewrite it as follows: 

Z^EUl^+^i^^.^.+^nr,)] (3.4) 

with V = A/(l — A). The random cluster representation is generated by eval- 
uating the product over all bonds in Eq. (3.4) and then summing over the 
possible values of a's and r's. If we think of the bonds originating from the a 
variables as green (g), and the bonds originating from the r variables as orange 
(o), then each of the resulting terms in the partition function can be graphi- 
cally represented as a collection of green and orange clusters as well as empty 
sites. The clusters of different colours may share sites, but not bonds. Denot- 
ing the configurations of green and orange bonds as ujg and cJq respectively, 
we can then write the partition function as 

Z (xJ2 t;^(^g)t;^(^°)2^('^^)2^('^°) (3.5) 

with b{u) being the total number of bonds (of a specified colour), and c{uj) 
being the number of corresponding connected components. The rule for count- 
ing the connected components is as follows: every site that is not a part of a 
green cluster is considered to be a separate connected component for the pur- 
poses of c(cjg), even if this site is a part of an orange cluster, and vise versa. 
In particular, the quantities i;''('^g)2'^('^s) ^^^1 y^('^°)2'^('^°) are to be interpreted 
exactly as in the usual random cluster models. 



3.3 Self-duality and criticality at X = 1 



The duality relations for such random cluster representation of the standard 
AT model in two dimensions were established in Refs. [22] and [23]. Firstly, 



starting from any plaquette - a particular choice simply determines the overall ad- 
ditive constant. In this sense the mapping between the present random surface and 
the loop model is nonlocal. 
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let us write the generic AT Hamiltonian as 



-(3H = Y.[k (6,^,^ + 5.,.J + L6,,.^6^^,] . (3.6) 



The graphical representation for the partition function is then obtained along 
the lines of the previous section. The only difference is that this time double- 
coloured (i.e. green and orange at the same time) bonds are also allowed. The 
graphical weight of a given bond configuration oo is 

Wiuj) = /[K^g'^^o) j^b{uigAuio)24^g)2^^'^°^ (3.7) 

where 



eL+2K _ 2e^ + 1 
e^-1 



A = e'^ -I and B = ^ . (3.^ 



Observe that Eq. (3.5) describes a particular case of this model provided that 
A = v = A/(l - A) while 5 = 0. 

The dual model is obtained by placing orange bonds between the sites of a 
dual lattice every time when it does not cross a green bond of the original 
lattice. Correspondingly, the green bonds on the dual lattice are dual to the 
original orange bonds. The duality relations are given by: 

A* = 2B-^ and B* = 2A~\ (3.9) 

And the model becomes self-dual when AB = 2. It is therefore suggestive that 
our model becomes self-dual at A = 1 (or v = oo). In order to show that it is 
indeed exactly self-dual, we shall perform the above duality transformation to 
the orange bonds only (cUq -^ ^o); leaving the green bonds intact. This results 
in having green bonds on both original and the dual lattices. The green bonds 
on the dual lattice can be then split into those traversal to the original green 
bonds and those traversal to the previously vacant bonds, or symbolically 
u* = Qg y Q^. Here cu* is the configuration of (green) bonds dual to the 
orange bonds while fig and Q0 are the configurations of bonds transversal 
to the original green and vacant bonds correspondingly. The corresponding 
weight is now given by: 

We now observe that b{Qg) = b{ujg), and also that ^(^0) ^ as f ^ 00 (the 
original random cluster model becomes fully packed according to Eq. (3.5)). 
Then the weight in this limit becomes simply 

^r _ 2^('^g) 2'^(^g) 2'^''^*^-' (3 11) 
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The model described by such weights is manifestly self-dual since u;* = fio V 

Q0 -^ Qo and Q* = UoM oo^ -^ cUq as f — > 00, and there is a symmetry between 
the green and the orange bonds. 

It is tempting to speculate that a phase transition occurs exactly at the self- 
dual point. Although this is plausible, it is not the only possibility. In par- 
ticular, there may be a phase transition at some At < 1. However, we can 
say the following: If at A = 1 there is no magnetisation (i.e. percolation of 
green or orange bonds) then the theorem proved by two of us [24] applies; 
A = 1 is a critical point in the sense of infinite correlation length and infinite 
susceptibility. The only other possibility is positive magnetisation at A = 1 
which implies a magnetic transition - which could be continuous or first order 
- at some At < 1. (In particular, this is shown to happen, with a first order 
transition for the large-g versions of these models [25]). Although we find these 
alternative scenarios unlikely, we have, in any case, established the existence 
of a transition in this model for some value of A between zero and one. 



3.4 Speculative remarks on relation to the critical 4-state Potts model 



As mentioned above, in our opinion the most likely scenario is that a phase 
transition occurs precisely at A = 1. The interesting question then is that of 
the universality class. Without any supporting mathematical statements, we 
suggest that at this point our model behaves similarly to the 4-state Potts 
ferromagnet at its critical point. In order to substantiate this claim, let us 
first recall the random cluster representation for a g-state Potts model: 

Z = Y,K^^'^\<^\ (3.12) 

For g < 4 this model is universally accepted to have a continuous transition 
at the self-dual point K = ^g. Thus for the g = 4 model at the self-dual point 
we have 

Z = ^2^('^)4'=('^). (3.13) 

UJ 

On the other hand, we can use Eq. (3.11) to rewrite the partition function of 
our model at A = 1 as follows: 

Z ocJ2 2^^""^ 4^^"") 2^(^)-^(^). (3.14) 

The difference between the two models is in the last factor of 2'^(^)^^('^) in 
Eq. (3.14). It is, however, reasonable to speculate that it can be neglected. 
Indeed, on average c(f2) = c{u), and therefore one would expect the typical 
value of the difference c{Q) — c{uj) to be sublinear in the system size. By 
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contrast, the individual terms c{Q) and c{uj) indeed scale linearly with the 
size of the system so this correction may be "unimportant" . 

This, however, does not mean that the two models approach the self-dual point 
in a similar fashion. In other words, we expect the exponents associated with 
the critical point itself (such as rj and 6) of the two models to be the same, 
while this needs not be true for the exponents associated with the approach to 
the critical point (such as a, {3 and z/). In fact, the A = 1 point of our model 
may well be an edge of a critical Kosterlitz-Thouless phase in which case the 
approach exponents would take on extreme values (zero or infinity). 



4 Reflection positivity and phase transitions in the large n limit 

4-1 Reflection positivity 



This section concerns the reflection positivity property of the loop models 
defined by Eq. (1.1) which in turn permits the analysis of their large n limit. 
Let A denote a rf-dimensional torus. Here, and for the remainder of this paper, 
it will be assumed that the linear dimensions of A are all the same, and are 
of the form L = 2^. We denote by A^ = L'^ the number of sites in the torus. 
Let Q denote the set of all possible loop configurations on A. Finally let V 
denote a hyperplane perpendicular to one of the coordinate axes which cuts 
through the bonds parallel to this axis dividing the torus into two equal parts. 
Let Qc and Qti be the bond configurations on the two sides of the "cut" , with 
the bonds intersected by V belonging to both sets. Thus Q = Q^^U Q-n, while 
Qp = Q^ n Qti contains only the intersected bonds. We now define a map 
'&V '■ Gc ^^ Gn such that it simply reflects the configuration on the "left" to 
that on the "right". Let / : Q-ji ^ M be a function that depends only on 
the bond configuration on the right and define '&pf{Qc) = /(''^(fl-c)) for any 
dc £ Gc- Similarly, we can use i^-p to map Q-ji -^ Qc, with this convention 'dp 
is the identity. A probability measure /i on the set Q is said to be reflection 
positive if for every such V and any functions / and h as described above 
{f^vf)^. > and {h ^f), = (/ ^ph)^. 

Theorem 4.1 The measures n determined by the weights in Eq. (1.1) are 
reflection positive on any even d-dimensional torus. 



Proof. Let A denote one such torus and V denote one of the above described 
planes. Let Qp E Qp he a configuration of bonds going through this plane. 
Assuming that fii^g-p) ^ 0, let us consider the measure /i(- | 0p). Our claim is 
that this splits into two measures, which we will call //£(■ | %p) and /i7j(- | 0p), 
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defined on Qc and Q-ji which are independent and identical under the reflection 

Indeed, it is not hard to see that /i(0p) 7^ if and only if Q-p has an even 
number of bonds of each colour. In each half of the torus, the endpoints of 
these bonds serve as "source /sinks" for bond conflgurations. In other words, 
a conflguration in, say, Qc must contain lines of the appropriate colour that 
pair up these sources. But, aside from having to satisfy these "boundary con- 
ditions", the weights are the same as given in Eq. (1.1). These two measures 
deflned accordingly on Qc, and Q-ji are the above mentioned ^c{' I Qv) ^"^^ 
fJ'Tzi- I dv) respectively. 

It is clear that if Qc ^ Qc then 

fJ'C {Qc I 9v) = l^n {^{Qc) I Qv) ■ (4.1) 

Furthermore if g-p is a conflguration and Qc is any conflguration that agrees 
with Q-p and has non-zero weight then for every 07^ G Qtz we see that 

l^niBn I 9v) = Kdfi I Qc) ■ (4.2) 

Thence, for every / that is determined on Q-ji we have 



SV 
8V 



(4.3) 
which cannot be negative. Similarly we get that {h 'dpf')^ = (/ 'dph) ^. □ 



One of the important consequences of reflection positivity is a Cauchy-Schwartz- 
type inequality: 

(/ ^vh), < ^{f^f),{h^h), . (4.4) 

which in turn leads to the chessboard estimates to be described below. (The 
reader interested in a more detailed description of reflection positivity is re- 
ferred to the review [26] and the references therein). 



4.2 Uniform exponential decay for large n 



In this subsection we will consider some n-colour models with n ^ 1 and 
vertex factors that are uniformly bounded above and below independently of 
n: < c < ui, . . . Um < C . Examples include the 0(n)-type models and the 
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corner-cubic models discussed in Section 2. However, the face-cubic model does 
not fall into this category since all the multi-coloured vertex factors vanish. It 
is no coincidence that we cannot treat these models since, as is obvious such 
models have a colour-symmetry broken phase for high enough value of bond 
fugacity (for brevity we omit a formal proof). 

The suppression of long contours will be established by showing that long lines 
of any particular colour are exponentially rare in the length of the line. To 
prove this we will need the so called chessboard estimate which in the present 
context reads as follows: 

Proposition 4.2 For x G A let cui, . . .oj^ denote indicator functions for bond 
events that are determined by the bonds emanating from the site x. (The Uj 
need not all be distinct.) For any of these ujj, cover the lattice with (multiple) 
reflections of the corresponding event and let Zj denote the partition function 
constrained so that at each site, the appropriately reflected event is satisfied. 
Then 

(n-,w>.<f #)*...(#)*. (4.5) 



Z J \z 



Proof. See Section 2.4 of Ref. [26]. 



Our principal result of this subsection: 

Theorem 4.3 Consider an n-colour loop model as described by Eq. (1.1) on 
the torus A (which is taken to be "sufficiently large") and suppose that the 
vertex factors are bounded below by c > and above by C < oo uniformly in 
n. For sites x and y in A, let Cx,y denote the probability that these sites belong 
to the same loop. Then, provided n is sufficiently large, there is a ^n > ^ such 
that for all values of R, 

C^,y < i^e-l^-^l/^" (4.6) 

where \x — y\ denotes the minimum length of a walk between x and y and K 
is a constant. 

Remark. For conceptual clarity, we will start with a proof of the case d = 2; 
all of the essential ideas are contained in this case. The problems in d > 2 
involve some minor technicalities and the general proof can be omitted on a 
preliminary reading. 



Proof of Theorem 4.3 {d = 2). Let us focus on a particular colour - red 
- and show the statement is true for red loops; this only amounts to a factor 
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of n in the prefactor. We define tlie "red event" as an event where at least 
two red bonds are connected to the site in question. It is clear that there are 
two main types of red events: those where the red bonds attached to a given 
site line up along the straight line and those where they form the right angle. 
These two types are shown in Figures 4(I-a) and 4(II-a) respectively. We will 



a 



II 



n: 



m 



r-r 



Fig. 4. Chess-board estimate on the "red" events of types I and II. The events 
are shaded in (a). Parts (b) and (c) represent the results of the first two reflections 
with respect to the dashed lines. The resulting tilings of the entire plane (torus) are 
shown in part (d). 



denote by uj 



(a) 



a = 1,2 the events of the first type and uJjj , (3 = 1, ... 4 those 
of the second type. Obviously there are only two distinct constrained partition 
functions which we respectively denote by Zj and Zjj. Figure 4 shows these 
single events (a), the results of the first two reflections with respect to the 
dashed lines (b), (c) and finally the configurations obtained by applying each 
reflection id = log2 N times in order to completely tile the surface of the torus. 
The grey lines correspond to yet unidentified bonds - these are the degrees of 
freedom left after the process of tiling has been completed. 

Let us perform the estimate on Zjj/Z first. We claim that if t// is any legiti- 
mate configuration that contributes to Zjj each of the red squares - of which 
there are N/4 - can be independently replaced by vacant bonds or a square 
loop of any other colour or left as red. Of course this may cost us an exchange 
of the "best" for the "worst" vertex factor but even so, the result is 



'II 



< 



R^ 



l + ni?4 



C 



N 



< 






N 



(4.7) 



Let us now turn attention to the Zi estimate. We start with a factor of R^ for 
the red bonds already in place - as well as another worst case scenario of C^ . 
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As for the lines that are orthogonal to (horizontal in Fig. 4(I-d)) once started 
in any colour they must continue until they wrap the torus, a total length of 
L = yN . There are n possible choices of colour for each line as well as the 
possibility of no bonds at all. Since there are a total of L lines altogether, this 
gives 

Zj<C^i?^(l+ni?^)^. (4.8) 

To obtain our estimate on Z, we simply pick the even (or odd) sublattice of 
dual sites and surround each site with one of n coloured elementary loops 
or with a "loop" of vacant bonds. Folding in the worst case scenario for the 
vertex factors this gives 

Z>c^(l + ni?^)T. (4.9) 

The ratio may be expressed as a product of two terms namely [R/ {l+nR^y/'^]^ 
and (1 + nR^)^/{l + ni?^)^/^ - times an additional {C/c)^ . Clearly the first 
term is bounded by n"^/^. As for the second ratio, if i? < 1 we may neglect 
nR^ for N ^ 1 and the ratio is bounded by one. On the other hand, if i? > 1, 
we may neglect the 1 and we get, modulo a factor of n^ ^ another n~^^*] we 
will settle for the bound of 1. 



Thus we have 



Z,\T! ( \ C 



lim f "< A^ (4.10) 

with the same upper bound for (^)^^^ valid for all A^. We thus denote the 
mutual upper bound by e„ ~ •nT^I'^. The desired result now follows from a 
standard Peierls argument: If x and y are part of the same loop, some subset of 
this loop must be a self- avoiding walk of length at least 2\x — y\. We enumerate 
all such walks and use the chess board estimate on each particular walk. Then 
if e„A2 < 1 where A2 is the two dimensional connectivity constant we write 

e„A2 = exp{-l/2U (4.11) 

(the factor of 2 because we must go there and back) and the stated result 
follows. □ 



Proof of Theorem 4.3 (rf > 2). The preliminary steps are the same as the 
two-dimensional case: There are again only two types of o^'s (but with more 
indices) and two constrained partition functions which we again denote by Zi 
and Zii. The pattern for Zn is the two-dimensional pattern in Fig. 4(II-d) 
reflected in all directions orthogonal to the plane visualised. Noting that each 
reflection doubles the number of red squares, we see that in the Zn patterns 



17 



there are a total of (2^)^^"^ x jL^ = jN squares altogether. Repeating the 
argument leading to Eq. (4.7) we end up with exactly the same bound. What 
is a little harder is the estimates on Zj and the partition function itself. 

We first claim that Z can be estimated by 

Z>c^(l+ni?^)^^. (4.12) 



To achieve this, we assert that the following holds: There is a set of plaquettes 
of the lattice with the property that each bond of the lattice belongs to exactly 
one plaquette. Once this claim is established it is clear that Eq. (4.12) holds; 
indeed there are just [1/4] x Nd plaquettes in question, we consider those 
configurations in which each of them is independently left vacant or traversed 
with an elementary loop in one of the n possible colours. Let us turn to a 
proof of the above assertion. 

Let §1, . . . e^ denote the elementary unit vectors. We adapt the following nota- 
tion for plaquettes: If, starting at x G Z'^ we first move in the e^ direction then 
in the e^ direction and then complete the circuit we will denote this plaquette 
by [ejOej.]x In general we can have [±ejO±efc]x and it is noted that [ejOe^Jx 
= [efcOe^Jx 

Starting with the origin, consider the following list of instructions for plaque- 
ttes: 

[eiO-e2]o, [e2 0-e3]o,...,[erfO-ei]o. (4.13) 

So far so good - each bond emanating from the origin belongs to exactly one 
plaquette. If x = {xi,...Xd) let us define 7rj(x) = (— 1)^^ to be the parity 
of the j coordinate. Then at the site x, select the following: [7ri(x)ei o 
— 7r2(x)e2]x, • • • [7rrf(x)erf o — 7ri(x)ei]x. Taking the union of all these lists, it is 
clear that indeed each bond belongs to at least one plaquette, the only question 
is whether there has been an over-counting. 

To settle this issue we establish the following assertion: Let x denote a site. 
Then a plaquette is specified by the instructions at x if and only if the same 
plaquette is specified by the neighbours of x on that plaquette. 

Indeed, consider one such plaquette namely [7rj(x)ej o — 7rj+i(x)ej+i]x. (If 
necessary we use the convention e^+i = §1). One of the relevant neigh- 
bours is x' = X + 7rj(x)ej. But at x', we have instructions for the plaquette 
[Kj{x')ej o — 7rj+i(x')ej+i]x'. Since vrj(x') = — 7rj(x') and 7rj_|_i(x') = 7rj+i(x') 
this is seen to be the same plaquette. The other neighbour, at x = x— vrj+iej+i, 
follows from the same argument. Reversing the roles of x and x' (as well as x 
and x) establishes the 'only if part of the assertion. 
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It is thus evident that no bond can belong to two plaquettes: If x and x' are 
neighbours and there are instructions coming from x to make some particular 
plaquette that includes the bond (x, x') then: 

(a) These are the only instructions coming from x pertaining to this bond. 

(b) The list at x' (and the other corners of this plaquette) also include this 
plaquette. 

(c) If y is another neighbour of x which is not in this plaquette there cannot 
be instructions at y to include the plaquette containing x, x' and y. 

Items (a) and (b) are immediate. To see item (c), note that by the assertion, 
such instructions would also have to appear on the list at x. But, by item (a), 
they don't. Equation (4.12) is now established. 

Finally, let us estimate Zj from above. First it is noted that the Zj patterns 
consist of parallel lines (pointing in the direction of the original red bonds) 
that pass through every site. Hence the only loops we can draw are confined to 
the various (rf — l)-dimensional hyperplanes orthogonal to these lines. Modulo 
vertex factors - which we estimate by C^ - each hyperplane yields the (rf — 1)- 
dimensional partition function. We thus have 

Zn < i?'^C^(H)^ (4.14) 

where S is the {d — l)-dimensional partition function with all vertex factors 
equal to unity. (This observation is not of crucial importance - the key ingre- 
dient is the number of available bonds - but it serves to compartmentalise the 
argument.) Let us tend to the estimate of S. 

We denote by Sm the contribution to S that comes about when there are 
exactly M bonds in the configuration. Given the placement of the bonds, 
let us count (estimate) the number of ways that they can be organised into 
self-returning walks and then independently colour each walk. Throwing in a 
combinatoric factor for the placement of the M bonds and the fact that there 
cannot conceivably be more than M/4 loops to colour we arrive at 



^n- X Y{M) (4.15) 



where Y{M) is defined as follows: For M bonds, consider their placement on 
the lattice chosen so as to maximise the number of ways in which the resulting 
(colourless) graph can be decomposed into distinct loops. (Such decomposi- 
tions are done by "resolving" all intersecting vertices into connected pairs of 
bonds, similar to Fig. 1.) Then the quantity Y{M) denotes this maximum 
possible number of decompositions. 
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Let 3 (which equals 2{d — l) — l = 2d — 3) denote the greatest possible number 
of local options for an on-going self-avoiding walk. We claim that 

Y{M) < (2^)^. (4.16) 

Indeed, take this optimal placement of M bonds and, starting at some prede- 
termined bond (and moving in some predetermined direction) draw a loop of 
length P. There are no more than 3^ possibilities for this loop. Then what is 
left over cannot yield a total better than Y{M — P). This must be done for 
every possible value of P and summed. Defining Y{0) = 1, we arrive at the 
recursive inequality 

Y{M)< J2 d^Y{M-P) (4.17) 

0<P<M 

and the bound in Eq. (4.16) can be established inductively. Putting together 
equations (4.15) and (4.16) we find 

- < E [^'^~M^'~')R''n^{2^r = {l + 2^Rn^r''^^'-\ (4.18) 

And thus 

Zi < C^i?^[(l + 23i?ni)('^'^)^''Y = [CR{1 + 23i?ni)(^-i)]^. (4.19) 



An estimate for the straight segments is nearly complete. We write 

.Z\ - C (1 + nR^I^ • ^ ' 

Let us split the power of the denominator: d/A = 1/4 + {d — l)/4; the first 
term will handle the R in the numerator and this ratio is less than n~^/^. As 
for what is left over, it is easy to see that 

1 + 2i,Rn^ 1 Rn^ 

—^^ r = r + 23 r<l + 23. 4.21 

(l + ni?4)i (l + ni?4)3 ^(l + ni?4)3- 

Thus we have 



1 
Z 



'I 



"< -(1 + 23)^^-^41 (4-22) 



c n 



1/4 



so again all of the cj's have estimates with the scaling of n~^^'^. The remainder 
of the argument follows the same course as the two-dimensional argument with 
modifications where appropriate. □ 
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4-3 Translational symmetry breaking and related phase transitions 



So far we have shown that for a sufficiently large n, there are no phase transi- 
tions associated with the divergent loop correlation length. However, this does 
not rule out the possibility of phase transitions of a totally different nature 
- in fact we shall see that such transitions do indeed take place. The type 
of these transitions happens to be sensitive to the lattice structure and the 
vertex factors that correspond to loop intersections. This leads to a variety of 
critical phenomena which has not been observed in the case of a honeycomb 
lattice [7]. 

We shall consider the case of all vertex factors being positive and uniformly 
bounded. Let the number of colours n be very large and consider the high 
fugacity limit of i? ^ oo. The model becomes fully-packed in this limit, i.e. 
every bond is occupied. But since n is large, according to the estimates of the 
previous subsection the partition function is dominated by configurations in 
which there are as many loops as possible. In two dimensions, it is clear that 
these "maximum entropy states" break the translational symmetry: the loops 
of different colours may run around either odd or even plaquettes of the lattice. 
(We remark that in rf > 2, the situation is considerably more complicated and 
is currently under study.) 

On the other hand, if the bond fugacity is low, the system must be transla- 
tionally invariant. Therefore one would anticipate that at some large but finite 
value of i? a translational symmetry-breaking transition takes place. Due to a 
doubly positionally degenerate nature of the resulting high-density state, we 
find it natural to expect that such transition is of the Ising type. Indeed, it 
appears to be very similar to the transition in the Ising antiferromagnet. 

Let us now turn to the rigorous arguments supporting this qualitative picture. 

Theorem 4.4 Consider an n-colour loop model on 1? with vertex factors 
hounded below by c > and above by C < oo uniformly in n. Then, if n is 
sufficiently large there is a phase transition characterised by the breaking of 
translational symmetry. In particular, there exist e (sufficiently small) and A 
(sufficiently large) such that for R^n < e all states that emerge as limits of 
torus states are translation invariant. On the other hand, if R^n > A, there 
are (at least) two states with broken translational symmetry. 



Proof. Let us start with the formal proof of translation-invariance at low 
fugacity. 

Any given site belongs to zero, two or four bonds. Let 0, and denote the 
corresponding events and Zq, Zq and Z^ the constrained partition functions in 
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which every site is of the stated type. We will estimate Prob(0) < {Zq/Zq^ '" 






and Prob(©) < (Zq/Zq) '"; obviously Zq = 1. In calculating Zq we notice 
that in all configurations, each bond must be coloured and there are at most 
N/2 separate loops. This gives 



Z^ < 



N 



R^n^/^ Y{2N) (4.23) 



where Y{M) is the same quantity that appears in Eq. (4.15). Similarly, in Zq 
half the bonds are used (resulting in a factor of R'^) and there are no more 
than N/A loops. We get 



Zn < 



Rn"^ 



N 



Y{N) (4.24) 



essentially the square root of the above. 



In order to show that there is translation invariance it is sufficient to estab- 
lish the following: Let A and B denote local events, let Tx be the transla- 
tion operator to the point a; G A and let e denote any unit vector. Then we 
must show that for all x with |x| large the probabilities of /ia {A fl Ty^B) and 
/iA {A n Tx+e-B) are essentially the same. We note that either the supports of 
A and Tx-B are attached by a ^-connected path of sites of type Zq and Zq, or 
they are separated by a connected circuit of cites of type Zq. (Two sites are 
considered ^-connected if they share the same plaquette.) However, in case the 
support of A is surrounded by a connected circuit of empty sites, it is clear 
that for any e, A and T^A have the same probability. Thus we arrive at 

l/iA (A n Tx5) - /iA (A n Tx+e5) I < IIK {Ca,t^b) (4.25) 

where 'Ca,t^b is the event that there is a *-connected path of non-0 sites 
between the support of A and that of Tx-B. Now, provided Rii}/^ is small, 
A^A (Ca,Txb) ^ e~''l^l with e"** ~ aRn^^'^ for some constant a. Thence transla- 
tion invariance (among the states that emerge from the torus) is established. 

Let us now turn attention to the opposite limit, namely Rn^^^ ^ 1. We claim 
that with high probability the plaquettes on the lattice are of one of the two 
types: they are either surrounded by four bonds of the same colour, or by 
bonds of four different colours. 

We will again define site events representative of the two purported ground 
states. If a site is on the even sublattice, we say it is of the A-type if it has 
four bonds of two different colours, one colour occupying the positive ei and 
§2 directions and the other, the negative ei and e.2 directions. The odd site 
of the A-type is defined as the mirror refiection of an even A-type site. A site 
is said to be of B-type if the roles of the even and the odd sublattices are 
reversed - see Fig. 5. Any site which is not of one of these two types will 
be called "bad". There are three choices for such bad sites, namely housing 
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Even 



Odd 



A 



B 



Fig. 5. Even and odd sites of A- and B-types. 

no bonds, only two bonds and a "+" configuration in which the ei and — ei 
directions have the same colour and similarly for the ±§2 directions. The first 
two are the old and from the preceding portion of this proof, while the 
constrained partition function for the latter event, denoted (for the lack of 
better notations) by Z^ can be bounded by 



Z#<i? 



•2N^2y/N 



(4.26) 



On the other hand, we can always estimate Z from below by 
discussed previously. Thus we have (as |A|-^ 00) 

/iA(0)<i?-'n-l/2, 

/iA(e) < n-1/2 , 



Rn^l^ 



N as 

(4.27a) 
(4.27b) 
(4.27c) 



Thus almost all sites are either A-type or B-type. However, there is no con- 
straint that A and B sites cannot be neighbours. But, as we show below, this 
possibility is also suppressed for n ^ 1. Firstly, let us notice that the reflec- 
tions through bonds map the even into the odd sites (and vice versa). By 
definition, the A-even and A-odd events are mirror images, and similarly for 
B. Thus, under multiple reflections, A-sites map to A-sites and B's to B's. 
Hence, if (A-B) is the event that the origin is of the A-type and its right 
nearest neighbour is of the B-type, we have 



/^a(A-B) < 



fZp^. 



V z 



2/N 



(4.28) 



where .^a-b is the partition function constrained according to the pattern 
shown in Fig. 6. Now, the interior of the A or B columns still consists of the 
four-bond loops - just as a pure A or pure B-type tiling would. However, as is 
easily seen in Fig. 6, a boundary between the columns is formed by a zigzag 
line of a single colour running around the torus. Thus we arrive at 



>A-B 



R 



2N -N - 

n* 712 



^/n 



(4.29) 
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Fig. 6. Tiling of the torus obtained by multiple reflections of the (A-B) event. The 
sites marked as Ag are the even sites of the A-type etc. The zigzag lines that form 
loops wrapping the torus and separating the columns of A- and B-sites are shown 
black. 

and hence, in the large N limit, /iA(A-B) < n^^^"^. 

The rest of the proof follows easily: A connected cluster of A-type sites has, on 
its boundary a "bad" site or an A-B pair; similarly for the clusters of B-sites. 
Thus, if Rn^^^ and n are large enough, the possibility that any given site is 
isolated from the origin is small. Thence, given that the origin is of the A-type 
- which has probability very near one half - the population of B-type sites is 
small and vice versa. Evidently, there are two states, one with an abundance 
of A-sites and the other with an abundance of B's. It is not hard to see that 
in these two states, translation symmetry is broken - both are dominated by 
the appropriate staggered patterns. □ 



5 Final remarks and conclusions 



So far we have succeeded in proving the following general statements: 

• A generic multi-coloured loop model with all vertex factors i/i, . . . z/^ uni- 
formly bounded from above and below does not have a phase transition 
corresponding to the divergence of the loop size in any dimension provided 
that the number of colours n is sufficiently large. 

• In two dimensions these models undergo a different phase transition, pre- 
sumably of the Ising type, that is associated with breaking the translational 
symmetry^. While the examples of the models that have been considered 



^ At present it is not clear whether this statement holds in d > 2. While the "ground 
states" (the states at i? = oo) of the system are not translationally invariant, their 
degeneracy appears to grow very fast with d, and it is not obvious that the resulting 



24 



in this paper all have their vertex factors independent of the particular ar- 
rangements of different colours entering the vertex, all these results remain 
valid even if this is not the case as long as all different vertex factors are still 
reflect ion- symmetric and bounded from above and below by some positive 
numbers. 

A clear example of a model that does not follow this rule is the loop model de- 
rived from the face-cubic spin model in Section 2.3. Here the bonds of different 
colours are simply not allowed to share a vertex, making the corresponding 
vertex factor vanish. It is clear that as the bond fugacity R increases, the sys- 
tem will have a phase transition (possibly of the first order) associated with 
breaking the colour symmetry. Indeed, the only way to "pack" more bonds 
into the system is to force all of them to be of the same colour. This transition 
is very similar to the Widom-Rowlinson transition. 

Another interesting observation can be made about the two-dimensional loop 
model with all two- and four-leg vertices having the same weight (this is the 
model originating from the corner-cubic spin model - see Section 2.2). Com- 
paring the states of the system ai R = oo when all bonds are occupied (the 
fully packed limit), and at i? = 1 when no additional weight is associated with 
placing extra loops, we conclude that the {n + l)-coloured model at i? = oo is 
identical to the n-coloured model at i? = 1. Indeed, one should simply consider 
the vacant bonds at i? = 1 as being coloured grey to turn the n-coloured loop 
system into a fully-packed {n + l)-coloured system. 

Now start with the large enough n at R = cx). The result of Section 4.3 
guarantees the existence of a broken translational symmetry in this case. By 
the argument presented here this means the existence of a broken symmetry 
in {n — l)-coloured, R = 1 case. Now let us continuously increase the value 
of the bond fugacity in this (n — l)-coloured model. As R grows, there are 
only two possible scenarios: the broken translational symmetry is either lost 
via a phase transition or is retained all the way to i? = cxd. The former case 
corresponds to the intermediate symmetry-broken phase surrounded by the 
phase transitions at i? < 1 and R > 1. If the latter scenario is realised, then 
go to -R = oo and repeat the process of mapping it onto a lower n model. 
Notice that this latter scenario can not continue all the way to n = 1 because 
the n = 1, i? = 1 case is nothing but a loop representation for the Ising magnet 
at T = which does not have broken translational symmetry.[3 Therefore one 
must find an intermediate phase at some not very large value of n (although 
we cannot rule out a possibility of it being just the point R = 1). 



"entropy" will not destroy the transition at any finite value of R. 
^ Indeed, these loops appearing in the T = "high temperature" expansion are 
also the domain walls of the Ising model on the dual lattice at T — > oo. These dual 
spins are assigned their values independently with probability 1/2, and the contours 
separating regions of opposite type are manifestly translation invariant. 
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From the above we also learn that the n = 2 (Ashkin-Teller-like) model is not 
critical in its fully-packed limit (since it is just the T = Ising model). On 
the other hand the 0(2) loop model discussed in Section 2.1 is critical in this 
limit (it maps onto the square ice model). But the only difference between the 
two is the factor of 3 in the same-colour four-leg vertex factor. 

On the basis of this observation (along with the first two remarks of this sec- 
tion) we reiterate that the vertex factors and lattice details are often important 
for determining the phase diagram of a particular model. 



A Vertex factor for the 0(n) vertices and loop decomposition 



In this appendix we present the calculation of the generic vertex factor ap- 
pearing in the loop expansion of the 0(n)-type model featured in 2.1. As it 
turns out, this calculation is independent of the lattice details (such as the di- 
mensionality or coordination number) which is reflected in our notation. Thus 
let Qn denote the n-dimensional solid angle and let S = (5*1, . . . S'„) denote the 
components of a n-dimensional unit vector. Since we only consider the spin 
at one particular site, the actual site index is conveniently dropped for now. 
The vertex factor we need to calculate is 

^mi,...m„ = TTy-j- / 1 ^ . . . 0„ " di2„ (A.l) 

where |fi„|= 27r"'/^[r (n/2)]~^ is the total solid angle and dr2„/|r2„| is the Haar 
measure. The object ^'m.i,...m„ is, of course, the "vertex factor" originating from 
2mi terms from the first component, . . . 2m„ terms from the n-th component 
associated with these numbers of these types of bonds entering into a vertex. 

Proposition A.l Let (f){2m) - which is equal to {2m — 1)!! - denote the 
number of different ways of dividing 2m objects into m distinct pairs. Then 

i: mjT^O 

where M = X^iLi '^i ^^^ 

D„(M)=2-^r(^)/r(M+^) = (n-2)!!/(n + 2M-2)!!. (A.3) 

Remark. The calculation of the value for i'mi,...mn has appeared before, cf. 
Ref. [27] and references therein. However, the interpretation that emerges from 
this formula has, to our knowledge, not been discussed previously {cf. the 
remark following the proof). 



26 



Proof. The proof boils down to the calculation of the above integral. We 
write 

dfi„ = sin"-2 On d^„ dfi„_i (A.4) 

where 6'„ = S ■ e„. Also, 

S = (Tsin^„, cos^„) (A.5) 

where T is an (n — 1) -dimensional unit vector. Finally, we note the fact that 

0Fr(f-i) 
r2„ — r2„_i 7 — r — (A. 6) 

r(f) 



Starting with the definition Eq. (A.l) with the help of equations (A.4)-(A.6), 
the recursion relation is obtained: 






— f T^™! . . . T^^II^f df]„_i 

n—1 I -^iin— 1 

^y^' . f (sin ^„)'^- ™^+"-' (cos ^„)'™" d^n 



, 2 
X 



2 



2 



G^ 



where 






r(f 



Using the recursion relation A. 7 to reduce the dimensionality of the spin space 
to 1 (and the fact that u^i = 1), we arrive to the desired expression: 






2Mr(M + f) ,^^0 
n 



; X ^ ^ ^ (A.9) 

^'^ n (2m, -1)!!. 
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Remark. As a consequence of the expression (A. 2) we may decompose the 
coloured bond configurations to obtain a colourless loop model (which gener- 
alises the result of Eq. (2.3) obtained for the square lattice). In particular, at 
each vertex we have a factor which depends only on the total number of bonds 
entering the vertex (and the dimensionality of the spin n) times the number 
of ways that the bonds of the various colours can be paired. By choosing a 
particular pairing scheme at each vertex, the bond configuration breaks, un- 
ambiguously, into a collection of self-returning walks (loops). Configurations 
with the same set of loops - that is to say the same set of bonds and the same 
pairing schemes (walking instructions) at each vertex - but which differ in 
the colours of the loops are seen to have identical weights. Thus we consider 
configurations /C of loops (bonds + walking instructions) and the partition 
function becomes the analog of the expression 2.3: 




uZ'n'^'^^ (A.IO) 



with Up being the factor for a vertex with 2p bonds given by 



u„ = vJn) = , ^^ . = — ^ -, . (A.lll 

' ^^ ^ 2Pr(p+f) n{n + 2)...{n + 2p-2) ^ ' 



We observe that the right hand side of Eq. (A.IO) provides a well defined model 
for all n and can be continued - essentially with no ambiguity - to non-integer 
values. As usual, we may define correlation functions via loop probabilities or 
be content with the numerators containing strings divided by denominators 
without. 

In the previous discussions of these issues (see Ref. [27] and references therein) 
the goal has been to derive the self-avoiding walk as the n -^ limit of 
the 0(n) models. Thus the vertex weights themselves get continued without 
reference to intermediate models. (And indeed, the self-avoiding walk does 
emerge as n — *> 0.) Although we have only treated the case where the original 
spin models were defined hj Z = Trn(jj) (1 + ASj ■ Sj) as opposed to Z = 
Trn(jj) exp (ASj ■ Sj) - we claim that the latter case also leads to graphical 
models that are well defined for non-integer n. (The derivation is somewhat 
intricate and will eventually appear in a future publication.) However, we 
remark that in either case, the n — *> limit is, in this context, quite simple to 
understand. 



28 



Acknowledgements 



The authors would like to thank J. Kondev, R. Kotecky and T. Spencer for 
the interesting discussions. LC and KS were supported in part by NSA Grant 
No. MDA904-98-1-0518 and NFS Grant No. 99-71016; LPP was supported in 
part by DOE Grant No. DE-FG02-90ER40542. 



References 



[1] J. Kondev, Int. J. Mod. Phys. B 11 (1997) 153-9, |cond-mat/9607181 . 

[2] S.O. Warnaar and B. Nienhuis, J. Phys. A: Math. Gen. 26 (1993) 2301-16. 

[3] H.A. Kramers and G.H. Wannier, Phys. Rev. 60 (1941) 252-62. 

[4] M. Kac and J.C. Ward, Phys. Rev. 88 (1952) 1332-7. 

[5] N.V. Vdovichenko, Sov. Phys. JETP 20 (1965) 477-9. 

[6] R.J. Baxter, Exactly Solved Models in Statistical Mechanics (Academic Press, 
New York, 1982). 

[7] E. Domany et al.. Nuclear Physics 190 (1981) 279-87. 

[8] B. Nienhuis, Phys. Rev. Lett. 49 (1982) 1062-5. 

[9] B. Nienhuis, Int. J. Mod. Phys. B 4 (1990) 929-42. 



[lo; 

[11 

[12; 
[13; 

[14 

[is; 

[16 



[17 



[is; 

[19 



H.W.J. Blote and B. Nienhuis, J. Phys. A: Math. Gen. 22 (1989) 1415-38. 

M.T. Batchelor, B. Nienhuis and S.O. Warnaar, Phys. Rev. Lett. 62 (1989) 
2425-8. 

M. Aizenman and B. Nachtergaele, Commun. Math. Phys. 164 (1994) 17-63. 

M.J. Martins, B. Nienhuis and R. Rietman, Phys. Rev. Lett. 81 (1997) 504-7, 
|cond-mat/9709051 . 



R.J. Baxter, J. Phys. A: Math. Gen. 19 (1986) 2821-39. 

J. Kondev and C.L. Henley, Phys. Rev. B 52 (1995) 6628-39. 



J.L. Jacobsen and J. Kondev, Nuclear Physics B 532 (1998) 635-88, cond- 
mat/9804048[ 



J.L. Jacobsen and J. Kondev, J. Stat. Phys. 96 (1999) 21-48, cond- 



mat/9811085 . 



J.V. Jose et al., Phys. Rev. B 16 (1977) 1217-41. 
H.J.F. Knops, Phys. Rev. Lett. 39 (1977) 766-9. 



29 



[20] J. Ashkin and E. Teller, Phys. Rev. 64 (1943) 178-84. 

[21] D.B. Abraham and CM. Newman, Phys. Rev. Lett. 61 (1988) 1969-72. 

[22] L. Chayes and J. Machta, Physica A 239 (1997) 542-601. 



[23] C.E. Pfister and Y. Velenik, J. Stat. Phys. 88 (1997) 1295-331, |coiici 
I mat/97040T7[ 



[24] L. Chayes and K. Shtengel, Commun. Math. Phys. 204 (1999) 353-66, |cond 



mat/9811203 . 



[25] T. Baker and L. Chayes, J. Stat. Phys. 93 (1998) 1-15. 

[26] S.B. Shlosman, Russ. Math. Surv. 41:3 (1986) 83-134. 

[27] N. Madras and G. Slade, The Self-Avoiding Walk (Birkhauser, Boston, 1993). 



30 



